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Abstract 

Wc present a general formalism for the calculation of finite-width contributions 
to the differential production cross sections of unstable particles at hadron collid- 
ers. In this formalism, which employs an effective-theory description of unstable- 
particle production and decay, the matrix element computation is organized as 
a gauge- invariant expansion in powers of Fx/mx, with Fx and m,x the width 
and mass of the unstable particle. This framework allows for a systematic inclu- 
sion of off-shell and non-factorizablc effects whilst at the same time keeping the 
computational effort minimal compared to a full calculation in the complex-mass 
scheme. As a proof-of-concept example, we give results for an NLO calculation of 
top-antitop production in the qq partonic channel. As already found in a similar 
calculation of single-top production, the finite- width effects are small for the total 
cross section, as expected from the naive counting ~ Ft/mt ~ 1%. However, they 
can be sizeable, in excess of 10%, close to edges of certain kinematical distribu- 
tions. The dependence of the results on the mass renormalization scheme, and its 
implication for a precise extraction of the top-quark mass, is also discussed. 



1 Introduction 



The performance of the Large Hadron Colhder (LHC) so far has been extremely successful, 
with about 5fb~^ of integrated luminosity collected in the 2011 run and more than 20fb~^ 
in 2012. This has led to an unprecedented accuracy in the measurements of many Standard 
Model (SM) cross sections and distributions, both in the electroweak and strong sectors, 
and also to the discovery of a new particle consistent with a SM Higgs boson. With the high- 
energy run which will follow the 2013 shutdown, experimental errors are bound to reduce 
even further, leading to higher-precision measurements and to the possible observation 
of new-physics effects beyond the SM. This high experimental precision clearly motivates 
similarly accurate theoretical predictions for cross sections and kincmatical distributions, 
so that on the one hand the clean extraction of signals from the data is possible and, on 
the other hand, contributions to the backgrounds to processes of interest can be accurately 
constrained. 

Most of the phenomenologically interesting processes at the LHC, such as W and Z 
boson production, top-quark production and Higgs production, not to mention beyond- 
the-Standard-Model (BSM) processes, like supersymmetric (SUSY) particle production, 
involve massive unstable particles. These particles are not asymptotic states and show up in 
detectors as energetic jets and leptons, often being accompanied by large transverse missing 
energy ^t- The necessity of precise theoretical predictions therefore raises the question of 
how to correctly treat the decay of the intermediate unstable particle to the physical final 
states. For observables which are inclusive in the final states originating from the unstable- 
particle decay, it is often sufficient to treat the massive particles as stable, ignoring their 
decay. The error associated with this approximation is formally of order Tx/mx, where Fx 
and nix are the width and mass of the unstable particle. For the aforementioned processes 
this corresponds to less than a few percent, i.e. typically smaller than the experimental 
errors. Clearly, while the stable approximation is appropriate for the total cross section, it 
cannot be used to predict arbitrary kincmatical observables. Moreover, the error associated 
with this approximation could actually be numerically sizeable for new, yet undiscovered 
wide resonances, for example strongly-decaying SUSY particles. 

A step forward towards a realistic description of production and decay of an unstable 
particle X is the narrow- width approximation (NWA) which is a framework commonly used 
in the context of high-energy calculations for hadron colliders. In the NWA the particle 
is produced and allowed to decay to the physical final states while remaining on shell. At 
next-to-leading order (NLO), radiative corrections are given by factorizable virtual and real 
contributions to the on-shell production and decay subprocesses. While technically only 
slightly more involved than the stable-top approximation, the NWA preserves spin correla- 
tions between the production and decay subprocesses, and allows for realistic kincmatical 
cuts on the momenta of the physical final states (i.e. leptons and jets). 

The NWA includes neither off-shell effects related to the virtuality of the intermedi- 
ate unstable-particle propagator, nor non-factorizable corrections linking the production 
and decay subprocesses. Sub-resonant or non-resonant contributions, which correspond 
to diagrams with the correct physical final state but which involve fewer or no interme- 
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diate unstable-particle propagators, are also neglected. As in the stable approximation, 
these finite- width effects are expected to be small, of order Tx/frtX) for inclusive-enough 
observables. This is due to large cancellations between virtual and real non-factorizable 
corrections and also because of the suppression of non-resonant contributions. However, for 
arbitrary kinematical distributions and in particular, close to certain kinematical thresh- 
olds where the cancellations mentioned above are less effective, finite-width effects can be 
large. Strikingly, in Refs. [1-3] it was pointed out that the naive expectation of the error 
associated with the NWA can be underestimated by an order of magnitude for BSM pro- 
cesses where the mass of daughter particles approaches the mass of the parent particle X. 
This is relevant for searches of SUSY in decay cascades, where one often observes some 
degree of mass degeneracy between particles in different steps of the cascade. Recently, 
non-negligible off-shell effects were observed even in Higgs production and decay to massive 
vector bosons [4], due to interferences between resonant and non-resonant contributions. 
Thus, it is clear that an approach that goes beyond the NWA and which includes at least 
the dominant finite-width effects is desirable. 

A possible solution to the issue of finite-width effects is clearly the calculation of the 
full, gauge-invariant set of diagrams corresponding to a given physical final state. This 
approach includes the coherent sum of resonant and non-resonant contributions, treats 
the intermediate resonant particles as fully off-shell and contains both factorizable and 
non-factorizable corrections at NLO. Self-energy contributions can be resummed in the 
unstable-particle propagator in a consistent gauge-invariant way using, for example, the 
complex-mass scheme [5,6]. Examples applying the complex-mass scheme to production 
of unstable particles at NLO include the calculations of four-fermion production at an 
e+e" collider [5], Higgs decay to vector-boson pairs [7] and two recent independent calcula- 
tions of off-shell effects in ti production [8-10]. While the complex-mass scheme approach 
is completely general and very flexible, allowing the calculation of arbitrary kinematical 
distributions, the full NLO computation is technically challenging, requiring both the cal- 
culation of a much larger set of diagrams than for the corresponding on-shell process and 
special techniques to handle 5- or 6-point functions with complex masses. 

An alternative approach to the full NLO calculation was presented in Ref. [11] and 
applied to processes of t-channel and s-channel single-top production [11,12]. The approach 
of Ref. [11] is the generalization of the effective fleld theory (EFT) description of resonant- 
particle production of Ref. [13], which was employed in the calculation of inclusive H^-pair 
production at an e"'"e~ collider [14,15]. The EFT calculation results in a systematic, gauge- 
invariant expansion of the matrix elements in powers of Tx/^X) in a way which can be 
considered a generalization of the pole approximation [16, 17]. Compared to the full NLO 
calculation in the complex-mass scheme, the effective-theory approach has the advantage 
of identifying the terms that are relevant to achieving a given target accuracy prior to 
the actual calculation. This greatly reduces the complexity of the computation while at 
the same time allows for the inclusion of the leading off-shell and non-factorizable effects 
in a completely differential manner. For single-top production finite-width effects were 
found to be small for inclusive-enough observables, although they can be large close to the 
kinematical edges of some distributions. This general picture is consistent with the results 
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found by the full NLO calculations of top-pair production [18]. 

In this paper we give a second example of the application of the EFT formalism 
of Ref. [11] and calculate the cross section for the pair-production process gg — )■ — )■ 
W^W~bb. The top-pair production process has been studied extensively over the years in 
the stable-top and narrow-width approximations and, more recently, using the complex- 
mass scheme (a detailed list of references is given in Section 3) . It thus represents a perfect 
proof-of-concept calculation by which to test the validity of the EFT formalism, extend 
it to more than one unstable particle and to compare it to different available approaches. 
The paper is organized as follows: in Section 2 we review the effective-theory formalism 
and introduce a treatment of real corrections which differs slightly from the one used in 
Ref. [11]. The calculation of the LO and NLO relevant amplitudes for the specific exam- 
ple of ti production is described in Section 3. In Section 4 we present results for several 
distributions and assess the effect of finite-width contributions by comparing the effective- 
theory predictions with results obtained in the NWA. In that same section we also discuss 
the effects of using different mass-renormalization schemes (more precisely the pole and 
PS schemes) on the cross section and distributions, and their possible implications for a 
precise extraction of the top-quark mass from data. Finally, our conclusions are given in 
Section 5. 



2 Effective- theory description of unstable-particle pro- 
duction 

The effective-theory framework for the description of unstable-particle production used in 
this work was first formulated for the total cross section in Ref. [13] and applied to the case 
of inclusive W~^W~ production at an e~^e~ collider in Refs. [14,15]. The formalism was 
later extended to the more general case of differential cross sections and applied to single- 
top production at hadron colliders [11,12]. In this section we review the main features of 
this approach, referring the reader to the aforementioned references for further details. 

Unstable-particle effective theory is built upon the hierarchy of two scales, namely, the 
typical virtuality of the resonant unstable particle X, which is set by its decay width, 
pj^ — m\ ~ nix'^Xi and the particle mass mx- This hierarchy is encoded in the ratio 
'^x/^x ^ 1- The latter is treated as a small parameter, 5, in line with the strong and 
electroweak coupling constants a^, and and allows for a systematic expansion of the 
full matrix elements. These different expansion parameters are related by the counting 
scheme ^ 

2 /IN 

~ a ~ aeto . (1) 

mx 

In the following we will generically refer to any of the above parameters as 5. The expansion 
in 5 is implemented at the Lagrangian level, replacing the (B)SM fields and interactions 

^Note that in (1) we assume that the unstable-particle decay proceeds via electroweak decay channels, 
Tx oc aew For resonances decaying via strong interactions the counting scheme is Tx ~ as ~ \/oiew 
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with effective fields and vertices. The effective fields are associated with different mo- 
mentum regions defined according to the scaling of their momenta with respect to the 
parameter S and encode the physics at the two very different scales that characterize the 
production and decay process^. For the problem at hand these momentum regions are 
a hard region (go ~ |g| ~ fnx), a soft region (go ~ |g| ~ rnxS) and collinear regions 
{rii - q ^ mxS, Ui - q ^ mx, q± ~ mxy/S)- Here n,, -hi are light-like vectors associated with 
the momenta of the external massless particles, rii = (l,pj/|pj|), Ui = (1, —Pi/lPil), and q± 
is the remaining, perpendicular component of the momentum g. 

In the effective theory only low-virtuality modes with < are kept as dynamical 
degrees of freedom, and are described by effective fields in the Lagrangian. These include, 
in particular, a field $x to describe the resonant unstable particle X. The finite-width 
of the particle is resummed into the leading EFT kinetic term in a generalisation of the 
heavy-quark effective theory (HQET) Lagrangian in the case of a non- vanishing width [13], 

^^eItmu = '^^x^l [^^-9- ^) $x , (2) 

where rhxv, with f ^ = 1, represents the large, on-shell component of the resonant-particle 
momentum, and rhx is the renormalized mass in a generic renormalization scheme. The 
coefficient fix is related to the complex pole fix = 'm'x—i'mx^ x of the full unstable-particle 
propagator, 

Ux = /fi^^ . (3) 
mx 

In the pole scheme, rhx = f^x, has the simple form Qx = —"i^x- Additional terms 
in the EFT Lagrangian are given by bilinear terms for soft and collinear fields, power- 
suppressed corrections to (2) and terms describing the interaction of $x and collinear 
fields with soft fields. Hard modes are not explicitly part of the effective Lagrangian 
and their contribution is encoded in matching coefficients multiplying effective interaction 
vertices. These can be schematically parameterized as 

(4) 

where J-},, J-"^ and J^^ji denote functions of fields and derivatives and the indices i, j, k 
label different Lorentz structures. 0c,s generically represent collinear and soft fields and 
Ci^P and Cj^D are the hard matching coefficients of the production and decay effective 
vertices, which are computed from on-shell SM amplitudes. In this context "on-shell" has 
to be understood as Px = m'j^ + nix fix = A^x, meaning that the effective couplings in 
the Lagrangian are in general complex. This is a feature that the EFT framework shares 
with the complex-mass scheme. The interaction terms Ck^NR^xR encode the contribution 
of non-resonant configurations which also contribute to the cross section starting from a 
certain order in 5. Note that in order to describe pair-production of unstable particles. 



^Here and in the following we assume that the invariants Sij = 2pi ■ pj constructed from the external 
momenta are of the same order of rnj^, and we treat them as a single scale. 
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(a) (b) 

Figure 1: Feynman-diagram topologies contributing to the process iii2 — )■ /i/2- (a) reso- 
nant production through an intermediate unstable X; (b) non-resonant production. 

e.g. ti production, two distinct resonant fields, $4 and $f, have to be introduced which 
annihilate a top and an antitop state respectively. Furthermore, the kinetic terms will 
contain two velocities v and v which are generally different. 

2.1 Born amplitudes 

From a practical point of view, at tree-level the EFT result coincides with an expansion of 
the matrix element around the complex pole of the propagator. In fact, the EFT approach 
can be viewed as a generalization, to arbitrary order in S, of the pole approximation [16,17]. 
To be specific, let us consider the process 

H{pi)i2{P2) -> X{px), ... ^ h{ki)f2{k2) , (5) 

where the final-state particles /i and /2 can originate from the decay of an intermediate 
unstable scalar particle X (the generalization to higher spins is trivial), or from other 
production mechanisms. The Feynman-diagram topologies contributing to this process at 
the tree- level are given by the resonant and non- resonant contributions drawn in Figure 1, 
where the grey blobs denote the model-dependent production and decay vertices. The 
tree-level amplitude can be written as 

^^^^^ ^ VA{P.Ux)VnM}) ^ ^ 
Px - K,x 

where Vp and Vd represent the vertices for the production and decay of the (off-shell) 
particle X respectively, while M contains the non-resonant contributions. mo,x denotes 
the bare mass of the particle X. It is immediately clear from Eq. (6) that due to the inter- 
mediate resonant propagator, p\ — "^ox ~ ""^x^; only resonant diagrams will contribute to 
the amplitude at leading order in 6, while A/" will generally be suppressed by extra powers 
of S. The actual suppression of the non-resonant contributions is determined from the 
interplay between the suppression from the missing resonant propagators and the scaling 
of the couplings appearing in Vp, Vd and M, see Eq. (1). 
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The gauge-invariant expansion in S of the amphtude, Eq. (6), is obtained by expanding 
the matrix elements for production and decay around = fi\. This requires a projection 
of the external momenta {pi}, {ki} onto on-shell configurations {pi}, {ki}, with 

Pi = PiiPhPx) h = ki{ki,Px) ■ (7) 

The on-shell projection is chosen such that momentum is exactly conserved at each vertex 
and px = IJ'x- Note that the explicit form of the projection in Eq. (7) is not unique, being 
defined up to terms of order S. However, calculations at order 5" obtained with different 
projections, deviate from each other by sub-leading corrections of order 5""^^, i.e. always 
an order higher than the target accuracy of the calculation. 

Adopting Eq. (7), the expansion in S of the tree-level matrix element reads 



ApA. 



D 



tree 



A 



X 



where Ap = Vp({pi},px), Ad = Vd{Px-i {ki}) are evaluated with the projected momenta 
pi, ki, and the ellipses represent higher-order terms in 5. As expected, the leading-order 
amplitude (first line in Eq. (8)) corresponds to the resonant diagram with the vertices 
for production and decay of the unstable particle replaced by (gauge-invariant) on-shell 
amplitudes. These amplitudes are in fact directly related to the matching coefficients Co,p, 
Co,D of the (leading) effective vertices appearing in Eq. (4). The modified propagator 
Ax = p\ — ii\ = p\ — rh\ — rhx^x Dyson-resums the finite- width effects related to the 
self-energy Uxipx) of the particle X: 

11 

+ ... (9) 



p| - m| + Uxip^x) 

with higher-order terms in 5 indicated by the ellipses. Note that it is not the full self-energy 
that is kept in the propagator. Only the gauge-invariant hard part of the self-energy 
contributes to the matching coefficient and is resummed in the propagator. The sub- 
leading gauge-violating residual soft terms are included perturbatively and are combined 
with other contributions to form a separate gauge-invariant part of the one-loop amplitude. 
The next term in 5 in the expansion of the amplitude (second line in Eq. (8)) receives 
contributions from non-resonant diagrams, as expected, but also from resonant ones in 
which the propagator Ax is cancelled by terms of the form Pi — pi ^ S and ki — ki 5, 
which originate from the expansion around on-shell configurations. In the EFT language 
these contributions are described by effective four-particle operators, i.e. terms of the form 
Ck,NR^NR (4)- We stress that while the leading resonant term is gauge-invariant, 

at higher orders in 5 only the sum of resonant and non-resonant contributions is on the 
whole gauge independent. 
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Figure 2: Correspondence between the expansion by regions and the EFT calculation: 
the hard-region contribution (top left) corresponds to a 0{as) correction to the produc- 
tion matching coefficient Cj^p (top right), while the soft-region contribution (bottom left) 
reproduces one-loop soft-gluon corrections in the EFT (bottom right). 



2.2 Virtual corrections 

Virtual corrections in the EFT framework are divided into two categories: loop correc- 
tions to the EFT matrix elements, where the degrees of freedom flowing in the loops are 
given by the effective fields ^x, 0c, 0s, and corrections to the hard matching coefficients 
Ci,p, Cj^D, Ck,NR- The latter are computed by a "matching procedure" which arises by 
requiring that the full-theory matrix element and the EFT result coincide order by order 
in 5. They can be written as a power series in the couplings (the QCD coupling, in this 
case) , 



i,p 



c. 



(0) 

i.p 



Os 

27r 



(10) 



with similar expressions for Cj^o and Ck,NR- In practice, the EFT matrix elements and 
the matching coefficients can be obtained from an expansion of SM matrix elements using 
the method of regions [19], which we find to be more convenient here than an explicit 
two-step calculation of matching coefficients and effective-theory matrix elements. This is 
the approach we adopt for the calculation of the amplitudes for tt production in Section 3, 
where details on how the expansion by regions is implemented for the specific case of 
resonant top-pair production are given. 

As already outlined, the momentum regions relevant to the expansion are the hard, soft 
and collinear regions. The contributions of hard momenta (g ~ fax) correspond to the 
matching coefficients, while those from an expansion in the soft region {q ~ mx^) of the full 
SM integrals reproduces loop contributions in the effective theory. This is schematically 
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f 



f 



(a) 



(b) 



Figure 3: Examples of triangle and box diagrams contributing to the production and decay 
of the scalar X at one loop. 

depicted in Figure 2 for the case of the one-loop gluonic correction to the production 
vertex. Hard corrections correspond to factorizable contributions to the production or 
decay subprocesses [20], and, at leading order in S, coincide with the matrix element for 
production or decay of the on-shell massive particle(s). On the other hand, soft corrections 
encode non-factorizable interferences between the production and decay subprocesses as 
well as off-shell effects. Note that in general the contribution from collinear regions is 
needed to reproduce the full SM matrix element. However, for the case considered, this 
contribution vanishes if loop integrals are regularized dimensionally and the external light- 
fermion masses are set to zero, as done in this work. It can thus be safely ignored in 
the following discussion. In the soft region, the unstable-particle propagator in the loop, 
{px — qY — rrix ~ ^xS, is resonant and has to be Dyson-resummed, 



with the ellipses denoting, as usual, higher-order terms in 6. As in the tree-level matrix 
element, only the on-shell gauge-invariant part of the self energy is resummed in the reso- 
nant propagator. In the hard region, where {px — q)^ — ^'x ~ ^"x ^ ^x^x, no self-energy 
resummation is necessary and gauge-invariance is similarly preserved. 

An important feature of the effective-theory approach is that the EFT counting scheme, 
Eq. (1), enables one to assign the correct parametric scaling to any particular contribution 
of a Feynman diagram prior to its actual calculation. This determines whether or not it 
has to be computed to obtain a given accuracy in 6. 

Suppose for example that one wants to compute all terms which scale as as ~ 5^^"^ 
relative to the leading tree-level contribution. From now on we will refer to such corrections 
as NLO corrections. If we ignore the couplings appearing in the production and decay 
vertices Vp and Vd, the leading tree- level matrix element (first term in Eq. (8)) scales as 
1/Ax ~ and NLO is thus defined as the sum of all terms scaling as 5~^6^^'^ = 5"^/^. 
One can then immediately see that the sub- leading tree- level terms in Eq. (8), which 
scale as 6~^6 = 1 (assuming that VpVo ~ A^), are not part of the NLO matrix element, 
contributing to the amplitude only at NNLO. Consider now the triangle and box integrals 



1 1 



(11) 



(px - qy - rhj^ + '^xiiPx - q^) Ax-2px-q 



+ ..., 
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shown in Figure 3. In the soft region of the loop momentum (g ~ rrixS), the unstable- 
particle propagator inside the loop scales as while the gluon and fermion propagators 
scale as and respectively. Taking into account the scaling of the volume element 
d'^q ~ S^, one finds the scaling S~^as5'^5~'^S~^5~^ ~ for the soft part of the triangle, 
while the soft part of the box scales as as5^5^'^5^^5^^5^^ ~ 5^^/^. Therefore, the soft 
limit of both the triangle and box contribute at NLO. In the hard region (g ~ mx) all 
the propagators inside the loop scale parametrically as ~ 1 (in units of the mass mx)- 
Thus, the hard part of the triangle scales as 5~^asl^ /l^ ~ 5"^/^, while for the box one 
finds ttgl'^/l^ ~ 5^/^. In this case only the hard part of the triangle counts as NLO, while 
the hard box integral is highly suppressed, scaling as a N^LO correction. Given that for 
the simple process (5), the hard box in Figure 3(b) is the most complicated integral at 
one loop, discarding it simplifies the NLO calculation. As we will see in Section 3.3, for 
the phenomenologically relevant case of top-pair production one can use the same scaling 
arguments to show that all one-loop hard 5- and 6-point functions are parametrically 
suppressed and need not be computed to achieve NLO accuracy in the resonant region. 

2.3 Real corrections 

Even though the expansion by regions is well understood for loop diagrams, it is less 
clear how the expansion can in general be implemented for real corrections, since in the 
presence of an extra massless particle with momentum q, the expansion parameter can 
be p\ — jj?x, {px + l)"^ — fJ'X both. While the split into hard and soft contributions is 
possible for the total cross section, by relating real corrections to cut one-loop diagrams 
via the optical theorem, the formulation of the expansion in hard and soft contributions 
is not straightforward for an arbitrary observable. In Refs. [11,12] we proposed a way to 
circumvent this problem by using the exact real-radiation matrix element and expanding 
the integrated infrared and coUinear subtraction terms in 6 so as to properly match the 
singularity structure of the expanded virtual matrix element. 

While this "ad-hoc" treatment of the real matrix element was sufficient to illustrate 
the basic features of the EFT approach and to assess the contribution of non-factorizable 
corrections in single-top production, a theoretically more satisfactory treatment of real 
corrections is certainly desirable. In particular, one would wish for a framework in which a 
strict expansion in 6 for both virtual and real corrections is achieved and gauge-invariance 
is exact order-by-order in S. Here we introduce such a framework which will be used in 
Section 3.4 for the calculation of real corrections to tt production. 

Real gluon-emission contributions fall into three categories, as indicated in Figure 4. 
Topology (a) and (b) represent emission from the production and decay subprocess respec- 
tively, whereas topology (c) is due to generic non-resonant real-emission. Analogously to 
Eq. (8), the corresponding amplitude can be written as 

_ Ap+gAp ApAp+g 

"^real a ' \ ' "^sub-lead; \^^) 

where Axg = {px + qY — fJ'x ^P+g i^D+g) is the amplitude for on-shell production (de- 
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(a) 



(b) 



(c) 



Figure 4: Real gluonic corrections to the process iii2 — )• /i/2- 



cay) of the unstable particle with an additional real gluon. The first two terms in Eq. (12) 
contribute to the leading resonant real-emission amplitude, whereas w4sub-iead contains terms 
sub-leading in 5. Note that in this picture there is no topology with gluon emission from 
the unstable particle. In the corresponding Feynman diagram either of the two unstable- 
particle propagators can be resonant. In fact, when the emitted gluon is soft, q ~ mxS, 



both propagators can simultaneously be resonant, Px 



m, 



{Px + qf 



m, 



m\5. 



However, it can be shown [21] that the product of two unstable-particle poles can be written 
as the sum of two terms containing a single massive pole and a soft singularity 



(13) 



{p\ - m\x){{px + qf - ml^x) 2px ■ q 



Px - ^o,x (Px + qy - ml 



X 



so that the parameterization Eq. (12) still holds true. 

Following Eq. (12), the squared amplitude, A^reai = l-^reaiP, can be written as 



M 



real 



lAp+gAol'^ , \ApA 



+ ^ T"":'^ +2Re 



|A 



Xg\ 



{Ap+gAD){ApAD+g)* 



+ 



(14) 



The first two terms can be interpreted as factorizable real corrections to either production 
or decay of the unstable particle X, while the interference term gives rise to non-factorizable 
corrections. Note that the interference term gives a sizeable contribution only if both Ax 
and Axg are resonant, i.e. when the emitted gluon is soft. This shows that, for the real 
matrix element as well, non factorizable and off-shell effects are encoded by soft radiation. 
The omitted terms in Eq. (14), denoted by the ellipses, are suppressed by UgS ~ 5^/^ 
compared to the leading tree-level contribution and first contribute to the matrix element 
at N^LO in the power-counting. Hence, they can be neglected in most calculations relevant 
for hadron-collider phenomenology. However, we point out that they can, in principle, be 
computed to the desired accuracy in 5, in a way similar to the expansion of the Born 
amplitude, Eq. (8). 
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2.4 Cross sections 



The real and the virtual amplitudes squared have been split into a hard factorizable part 
and a soft non-factorizable part. The hard part is further divided into corrections to 
the production and corrections to the decay. As usual, the virtual and real cross sections, 
obtained by integrating the squared matrix elements over phase space, have infrared and/or 
collinear singularities. 

The cancellation of these singularities between real and virtual corrections takes place 
separately for all parts of the cross section. The (hard) production part typically has 
initial-state collinear singularities that have to be absorbed into the parton distribution 
functions. It can additionally develop infrared singularities and, when particles other than 
the unstable particle are present in the final state of the production subprocess, can develop 
final-state collinear singularities as well. These singularities cancel against the correspond- 
ing real singularities of the production part. The (hard) decay part typically has final state 
collinear singularities and infrared singularities that cancel against the corresponding real 
singularities. 

Finally, the non-factorizable soft corrections only have infrared, but no collinear singu- 
larities. Once again, these cancel between the real and virtual corrections. Also, with this 
cancellation, all explicit scale dependence vanishes and the non-factorizable corrections 
only depend implicitly on the renormalization and factorization scale. Because the soft 
emission is governed by a soft scale, /is ~ m^^, it is natural to choose a scale of this order 
for the coupling from the additional soft emission. The factorization scale and the scale 
related to the born term are not affected by this and should be kept at the usual hard 
scale. 

2.5 Mass scheme 

The results given in this section have been expressed in terms of the unstable-particle 
mass in a generic renormalization scheme, rhx, which is related to the gauge- invariant 
complex pole of the propagator by /i^ = + rhx^x- While the EFT approach is not 
limited to a particular scheme, it is more naturally formulated in a class of schemes related 
to the pole scheme by {rhx — fTix)/^x ~ ^, which have the property of preserving the 
EFT scaling p\ — m\ ~ m\5 under renormalization. The MS is not such a scheme 
since {rhj^^ — rnx)/rnyr^ ~ ~ ^A and the finite-width expansion becomes less 
transparent in this case. The pole mass on the other hand, which is the natural mass choice 
for unstable electroweak gauge bosons and other particles which do not couple to strong 
interactions, is affected by long-distance ambiguities related to QCD renormalons [22-24] in 
the case of heavy quarks. This makes a precise definition and extraction of this parameter 
problematic. It is possible to circumvent these issues with the use of so-called "threshold 
masses" [22,25-27] or "jet masses" [28,29], which preserve the EFT counting and are not 
affected by renormalon ambiguities. The implication of the mass-renormalization scheme 
choice in the case of ti production is discussed in Section 3.5. 
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3 Top-antitop production at hadron colliders 



At hadron colliders the main top-production mechanism is the production of a ti pair, which 
proceeds through strong interactions^. At Born level two partonic processes contribute to 
the cross section: qq — )■ ti, which dominates at Tevatron, and gg — )■ ti, which is the 
dominant production channel at the LHC Beyond leading order other partonic production 
channels contribute as well. Next-to-leading order QCD corrections to the tree-level results, 
both at the inclusive and differential level, have been known for more than 20 years [BO- 
SS] and have been recomputed recently for helicity-specific amplitudes using generalized- 
unitarity methods [S4]. One-loop electroweak contributions are also known [S5-S9]. At 
next-to-next-to-leading order in QCD several ingredients and partial results have been 
known for some time [40-51] and the full NNLO calculations for the fermion-initiated 
[52, 5S] and gq{q) [54] production channels have been recently completed. In addition, 
resummation of Coulomb singularities and soft logarithms up to next-to-next-to-leading 
logarithmic accuracy has been studied by several groups [55-59]. 

All the aforementioned results were computed in the approximation of stable tops. 
Results in the NWA, including the on-shell decay of the top-antitop pair to the physical 
final state W~^W~bb, were given in Refs. [60-6S]. The fully differential semi-leptonic decay 
of on-shell top quarks is also now known to NNLO [64,65]. Recently the full NLO QCD 
calculation of the off-shell top-pair production and decay process in the complex-mass 
scheme was performed by two separate groups [8-10]^. As expected, finite-width effects 
were found to be small for the total cross section. However, larger effects of up to tens of 
per cent were observed for more exclusive kinematical distributions [18], thus confirming 
the necessity of a careful treatment of finite-width effects both in ti production, and more 
generally, in processes involving unstable particles. 

In this section we present the calculation of the LO and NLO virtual and real amplitudes 
for ti production in the EFT framework, as a non-trivial example of application of the 
formalism described in Section 2. We will focus on the qq production channel, 

qq-^ti^ W^Whl 

where the decay of the two W bosons to leptons is understood. Our notation and con- 
ventions for momentum and colour labelling is given in Section S.l where we also give 
explicit expressions for the leading tree-level helicity amplitudes. The calculation of soft 
and hard virtual corrections as well as the renormalization necessary at NLO are detailed 
in Section S.S, while the implementation of real corrections is described in Section S.4. 

•^While the electroweak-mediated production of a single top is numerically only 3-4 times smaller than 
top-pair production, a much larger background makes the extraction of the signal in this channel more 
challenging. 

^The calculation of Ref. [9] treats also the leptonic W decays completely off-shell, while in Ref. [8] 
the W decays are described in narrow-width approximation. The effects of finite VF-width effects in ti 
production have been investigated recently in Ref. [10]. 
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3.1 Leading tree- level amplitude 

Throughout this work we adopt the following momentum labelling 

q{Pi; Ci)q{p2; C2) hhW^W~ h{pz\ Ci)h{pA\ C4)li{p5)MP(i)HP7)MP8) , (15) 

with C1...C4 the colour indices of the four external quarks. The leptonic decay of the two 
VF-bosons is treated in the NWA. For sake of simphcity, we introduce the shortcuts 

Pt = P3+P5+P6, 

Pi = P4+P7+P8, 

and the on-shell projections pt and pt, with pf = pf = fJ^f, which are necessary^ for the 
systematic expansion of the amplitudes in 6. As usual fi'f denotes the complex pole of the 
propagator, /ij = ml — irritTi. We also introduce the abbreviations 

/^t = pl-lA Ai = pl-fj,l (16) 

for the resummed inverse top and antitop propagator. The expression (16) is valid in the 
pole scheme, i.e rrit here denotes the pole mass. The conversion to other renormalization 
schemes is discussed in Section 3.5. The amplitudes for the tt production process are 
decomposed onto a basis of colour-state operators, given by 

o'^^^ - J_A A co^^^ - - (^7) 

^{c} ~ jy "C1C2'^C3C4 ^{c} — y]y2^7J C2C1-'C3C4 ' \^ ' J 

with T"^ denoting the generators of the SU (3) algebra in the fundamental representation 
and Nc = 3. The two operators describe the s-channel exchange of a colour- singlet and 
colour-octet state respectively, and satisfy the orthonormality condition 

Y^of^O^^S.,. (18) 

{c} 

At the Born level both doubly-resonant Feynman diagrams containing an intermediate 
top and antitop pair, as well as singly-resonant and non-resonant diagrams, contribute to 
the amplitude for the process (15). The dominant tree-level contribution to the amplitude 
is given by the leading term in the expansion in S of the left Feynman diagram in Figure 5. 
In the effective-theory framework, the interaction mediated by the intermediate s-channel 
gluon is replaced by a four-fermion contact interaction (right diagram in Figure 5) of the 
form 

• 2 

^Note that in an NLO calculation one can set pi = = TOj instead, since the term —initTt ~ mid is 
formally an NNLO contribution. 
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Figure 5: Leading tree-level doubly-resonant contribution to qq — ?> hhW^W^ in the SM 
and its EFT counterpart. The grey square and circles represent EFT production and decay 
matching coefficients. 

with s the partonic centre-of-mass energy, ai, q;2, a, a the Lorentz indices of the initial 
quark and antiquark and final top and antitop, and Ci, C2, c, c the corresponding colour 
indices. The tree level amplitude reads (omitting the on-shell decay of the two Ws) 

Aree = ^SZiog)^(O) , (20) 

with 

2 2 

= ^^v{P2)Y<Pi)u{p^W+{Pt + mt)^^{-pt + mt)f_v{P^) . (21) 

e± represent the polarization vectors of W^. Note that the numerator of A^^'^ depends on 
the projected on-shell momenta which guarantees the gauge-invariance of the result. The 
contribution of resonant and non-resonant diagrams will be discussed in Section 3.2. 
Using Eq. (1) one can easily determine the scaling of the leading amplitude, which is 

^(0) ^ ^ ^ r V2 . (22) 

As a consequence of the relative scaling of the top-quark width and the strong coupling 
constant, ag ~ ^JV^Jmt, the expansion of the amphtude is organized in powers of 5^/2. 
Thus in the following we define as N^LO all contributions to the amplitude that scale as 
(5{fe-i)/2 particular, NLO corrections, which we will focus on in this paper, scale as 
5°. Potentially, these include sub-leading tree-level terms, 0(a^ hard and soft virtual 
corrections to the leading doubly-resonant amplitude Eq. (20), and real QCD corrections. 
These three contributions will be discussed in turn in the following sections. 



3.2 Sub-leading tree-level contributions 

The first class of corrections which could potentially contribute at NLO in 8 is represented 
by terms arising from sub-leading tree-level diagrams, which are shown in Figure 6. These 
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(a) 





(c) 




Figure 6: Subleading tree-level diagrams in the effective theory: (a) singly-resonant dia- 
grams, (b) non-resonant diagrams, (c) doubly-resonant diagrams with electroweak match- 
ing coefficients (represented by the grey and white square in the figure). See the text for 
explanation. 

include QCD singly-resonant and non-resonant topologies (diagrams (a) and (b)). Since the 
effective-theory counting of Eq. (1) relates the expansion in as to the one in the electroweak 
coupling constant, electroweak diagrams must also be taken into account beyond LO in 6. 
For all practical purposes it is sufficient to consider the leading doubly-resonant topologies 
corresponding to diagram (c) in Figure 6. Here we use scaling arguments to determine the 
size of these terms compared to the leading-order amplitude, Eq. (21). 

Contributions with the structure of diagram (a) correspond to the terms on the second 
line of Eq. (8), where one of the resonant propagators is either cancelled by higher-order 
contributions originating from the expansion around the pole, or is not present in the first 
place. In the effective theory they are described by five-particle contact interactions, whose 
matching coefficients are proportional to g^gew Their contribution to the amplitude thus 
scales as 



(23) 



which makes them an NNLO correction. Their interference with the leading Born ampli- 
tude gives a NNLO correction to the cross section, while the square of diagram (a) con- 
tributes to the cross section initially at N^LO. Non-resonant QCD topologies correspond 
to six-particle effective interactions (diagram (b)) with matching coefficients starting at 
order glglw- Since they do not contain intermediate resonant propagators, these scale as 



^3/2 



(24) 



and contribute to the cross section as N^LO in 5, being suppressed by 5'^ compared to the 
leading Born amplitude. 

The doubly resonant electroweak diagrams (diagram (c)) scale as 



(25) 
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Figure 7: A sample of one-loop diagrams (above) contributing to the 0{as) corrections to 
the production matching coefficient (below). 

Note that in this case the suppression is given by the electroweak nature of the match- 
ing coefficients (~ g"^^) of the four-particle production operator, rather than by missing 
resonant propagators. The interference of diagram (c) with the LO amplitude thus gives 
an NLO correction. However, due to the colour structure of the electroweak production 
operators, this interference is non- vanishing only for a bb initial state. This contribution is 
numerically suppressed due to the smallness of the bottom-quark PDF inside the proton 
and is in practice negligible (< 0.1% of the leading doubly-resonant QCD-mcdiatcd dia- 
gram). The squared electroweak matrix element counts as an NNLO correction and in this 
case there is generally no suppression arising from the parton luminosities, since operators 
with any light-flavour quark in the initial state contribute. 

To conclude, if one neglects the (accidentally) tiny NLO contribution from the interfer- 
ence of diagram (c) with A^^\ the first Born-level corrections to Eq. (21) arise at NNLO 
from the interference of singly-resonant QCD diagrams with the leading amplitude and 
from the squared matrix element of electroweak doubly-resonant diagrams. Being sup- 
pressed by 5 ~ Tt/rrit compared to Eq. (21) they are expected to contribute to the cross 
section at the percent level. This is investigated more quantitatively in Section 4, in par- 
ticular to test whether the EFT counting applies in practice to the case of ti production. 
As we will see, sub-leading terms account for at most a few percent of the cross section 
which is in good agreement with the expectation from the EFT scaling arguments. 

3.3 Virtual QCD corrections 

As explained in Section 2.2, virtual corrections in the EFT approach are classified into 
corrections to the hard matching coefficients of the production and decay effective vertices 
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Figure 8: Two examples of one-loop diagrams which do not contribute to the hard matching 
coefficients at NLO in 6. See the text for further explanation. 

as well as soft corrections to the EFT matrix elements. The former encode physics at the 
large scale, q ~ nit, whilst the latter describe long-scale physics, q ~ mt6. To obtain NLO 
accuracy, corrections to the hard matching coefficients and one-loop soft corrections 
are required. As pointed out in Section 2.2, 0{as) hard matching coefficients can be 
related to one-loop corrections to the amplitudes for production and decay of on-shell 
tops. For the production part we use the results of Ref. [34], whereas the results for the 
decay can for example be found in Ref. [21]. A sample of the diagrams contributing to the 
production matching coefficient is shown in Figure 7. As for the toy model discussed in 
Section 2.2, hard loops connecting initial and final-state particles, or final-state particles 
originating from different decays, are suppressed by higher powers of 6. Consider, for 
example, the two diagrams shown in Figure 8. If the momentum flowing in the loop is hard, 
the intermediate antitop propagator in the first diagram is far off-shell. Thus the diagram 
is suppressed by a^Af ~ 5^/^ compared to the leading amplitude Eq. (21), which makes 
it a N^LO correction. The pentagon diagram, in which both top and antitop propagators 
are off shell, is even more suppressed, ~ Os^t^t ~ S^^"^, and contributes to the amplitude 
starting at N^LO. These contributions do not have to be included in a calculation that 
aims for NLO accuracy in 6, thus dramatically reducing the computational complexity of 
the EFT calculation compared to that in the complex-mass scheme. 

Contrary to the hard corrections, non-factorizable soft diagrams that connect initial and 
final states do contribute already at NLO in S, since the two intermediate top propagators 
remain resonant if the loop momentum is soft. The complete set of soft one-loop diagrams 
contributing to the amplitude at order 5° consists of: 

• 2 self-energy corrections, to the top and antitop resonant propagator, 

• 8 triangle diagrams, shown in Figure 9, 

• 6 box diagrams, shown in Figure 10, 

• 1 pentagon diagram, also shown in Figure 10 (right bottom corner). 

Note that the renormalization of the bare couplings and masses only affects the hard 
matching coefficients, but does not contribute to the soft corrections. The contribution of 
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the soft diagrams listed above can be obtained by expanding full QCD diagrams in the 
soft region. Only a consistent expansion in 6 guarantees the gauge invariance of the final 
result. Most of the integrals shown in Figures 9 and 10 were already computed for the 
case of single-top production [11, 12], with the new soft integrals required given by the 
triangle in Figure 9d, the boxes lOe and lOf, and the pentagon. As a further example of 
how the method of regions works we will now briefly discuss the calculation of the triangle 
diagram 9d. 

All soft integrals can be written in terms of the colour operators Eq. (17), the colour- 
stripped leading-order amplitude Eq. (21) and scalar functions depending on the invariants 
Sij = 2pi ■ pj built from the external momenta. In particular, the contribution to the 
amplitude of the triangle in Figure 9d can be written as 



3,d {a}- 

with the scalar integral lu defined as 



(26) 



[AnY J {2n)'^ At - 2pt ■ q At + 2pt ■ q 

and d = 4 — 2e. The simple expression of Itt follows from the eikonal form of the QCD 
fermion-gluon interaction vertex in the soft limit, u(p)7^(^f ±(^ + mf) ~ 2ptu{p), and from 
the expansion of the propagators in the integrand of the full QCD triangle integral, 

(pt ±qf - 111 = p1 ±2pf q + q^ - ^i] ^ At±2pf q, (28) 

where the quadratic term q^ ~ 5"^ must be dropped, being sub-leading compared to the 
term At ± 2pt ■ g ~ 5. Note that in Eq. (27) the off-shell top and antitop momenta pt and 
Pt have been replaced by the on-shell projections pt and pf, except in the combinations 
p1 — fif = At and pf — fi^ = At. In addition, finite-width effects have been resummed into 
the top propagators, in accordance with the EFT Lagrangian Eq. (2). 

Using Feynman parameters and standard integration techniques Eq. (27) can be reduced 
to the one-dimensional integral 

/«-=%e^''-r(l-e)r(2e) H dy [{y - e+)(y - ^r'^' f '^' ~ ^^^" "j '\ (29) 

27T Jq V ^tl^ J 

where we have introduced the invariant rj = 2pt ■ Pt/rrit and 

U = ^^V"''-'^'"- . (30) 
with io4. an infinitesimal and positive imaginary part. Expanding Itt in e one obtains 



(-1) 



1 / At 

- - 2 hi 



e \ TTitfi 



+ (31) 
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(a) 









Figure 9: Soft triangle diagrams contributing at NLO in the EFT counting. Dashed hues 
represent soft gluons. 



where the integrals I^^ and I^^"* are given by 



^ ln(-e+) - ln(-C- 



(32) 



AO) _ 



+ In(-e-) ln(e- - e+) + ^ ln^(^+ ' ' I " 



+Li2 



-21n(-e+)ln 



— Lh 



-In^ 



At + AfC^ 
At 



+ ln' 



At + AfC- 



. ^ , +21n(-^_)ln . ^ 



-2Li9 



At + Af^H 



+ 2Li2 



At-e- 



At + AfC- 



(33) 



Note that the only explicit scale dependence in Eq. (31) is due to the term In (— Aj/ {mtjj)), 
indicating that the natural scale choice for non-factorizable corrections is //soft ~ /S.t/mt ~ 
Ft. All other terms in Eq. (31) have no explicit // dependence and do not contain potentially 
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(a) 






Figure 10: Soft box and pentagon diagrams contributing at NLO in the EFT counting. 
Dashed hues represent soft gluons. 

large logarithms. This is a common feature of all soft contributions in Figures 9 and 10. 
The calculation of the missing box and pentagon integrals closely follows the one outlined 
above for the triangle and will not be given explicitly here. 



3.4 Real QCD corrections 

At NLO real QCD diagrams must also be taken into account. As discussed in Section 2.3 
and illustrated in Figure 11, any real emission diagram can be split into a linear combination 

of terms containing a single resonant propagator. In turn, these can be interpreted as 
corrections to either the production or decay subprocesses and their contribution to the 
amplitude is obtained from the expansion of the full QCD matrix elements around on-shell 
configurations. These on-shell configurations are defined by Pt = pj = IJ^I for diagram (a), 
{pt + ky — p^ — 111 for diagram (b) and pi = {pt + k)"^ = fif for diagram (c), with k being 
the momentum of the emitted gluon. Note that here only gluonic corrections to the tree- 
level topology are considered, but a full calculation which also includes the gluon-initiated 
production channel would require the calculation of diagrams with a gq{q) initial state at 
NLO. 

Soft and coUinear singularities are treated in the standard way, by adding a suitable 

subtraction term to the real matrix element and subtracting the corresponding integrated 
function from the virtual corrections, such that the cancellation of poles is manifest and real 
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Figure 11: Real-corrections diagrams contributing to the cross section at NLO in S 

corrections can be safely computed in four dimensions. In particular, for the calculation 
of the results given in Section 4, two independent numerical codes were used, one imple- 
menting the FKS subtraction method [66, 67] and the other, the Catani-Seymour dipole 
method [68,69]. As explained in Section 2.4, the cancellation of virtual poles against real 
ones occurs separately for factorizable and non-factorizable corrections. For factorizable 
corrections this follows from the cancellation of infrared and coUinear singularities for the 
on-shell production and decay subprocesses. On the other hand, the (infrared) poles of 
non-factorizable virtual corrections (i.e. the sum of all soft virtual contributions) are given 
by 
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(34) 



with Sij = 2pi-pj. These can be shown to coincide with the pole structure of the production- 
decay interference terms in Eq. (14). 



3.5 Mass-scheme choice 

The results shown in the previous sections have been computed in the so-called "pole-mass 
scheme", in which the renormalized mass parameter is chosen to coincide with the real 
part of the complex pole of the top-quark propagator. This particular mass scheme is 
affected by ambiguities related to QCD renormalons [22-24], making a determination of 
the top-quark mass at an accuracy better than Aqcd conceptually impossible. Avoiding 
these problems requires the use of a mass scheme which is not affected by long-distance 
ambiguities. An example of such a mass is the MS-renormalized mass, m^^. On the other 
hand the MS mass presents the unpleasant feature of not preserving the effective-theory 
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counting (or, in other words, it differs from the pole mass by terms of order as), thus 
spoiling the EFT counting scheme. More suitable choices are represented by the so-called 
"threshold masses" [22,25-27], introduced in the context of heavy-quark threshold physics, 
and "jet masses" [28,29] studied for boosted top-jet events in e"'"e~-collisions. These masses 
have the attractive feature of not being sensitive to low-momentum physics and of differing 
only by amounts of order ~ mjaf ~ rritS from the pole mass, thus preserving the effective- 
theory counting. In the following we will consider the "potential-subtracted" (PS) mass 
definition [25] as an alternative suitable mass scheme, illustrating how this fits into the 
effective-theory calculation. 

The relation between the pole mass and the PS mass can be found in Ref. [25], and up 
to second order in as (enough for applications in this work) reads 



m = mi,ps + /ips 



2n 



2n 



6. 



PS 



(35) 



with 



SI 



PS 



^PS 



2Cf 
Cf 



ai-&o(ln(^l-2 



(36) 



and Oi = 31/3 — lOnj/9, Bq = 11 — 2nf/3. The cutoff /xps has to be chosen of order ~ rritas 
to preserve the EFT counting. Substituting the pole mass in favour of the PS mass in the 
(renormalized) resummed top propagator one obtains 



(37) 



pf - ^i,ps + ^"^t,ps 



with the ellipses denoting terms that scale as 5^ or higher. Note that, counting /ips ~ 
the term — cts^P/ipsmi^ps/vr scales as the leading term pf —mfpg + init^ps^t ~ S, and has to 
be kept in the denominator of the propagator. The term proportional to af is suppressed 
by an additional \/5 and can be re-expanded to fixed-order in the expansion parameters. 



27r2 



Pt - mi,ps + ^^t,ps^t - v^i AiPS"^i,PS 1 Pt - mi,ps + imt,ps^t - v"^! AiPS"^t,PS 



.(3^ 



Note that the replacement rrit — )■ mt,ps + ^5f^/ips + ... in the numerator of the matrix 
elements does not generate extra terms at NLO since a^/ips ~ 5. Thus to convert the pole- 
scheme result for the renormalized cross section to the PS scheme it is sufficient to replace 
ml — irritTt — )■ ps — ii^t,ps^t + ^(^f^/^ps^^i.ps in the denominator of the top (antitop) 
propagator (in both the LO and NLO cross section), — )• ps in the numerator and to 
add the higher-order propagator contribution 



5X 
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polc-!>PS 



J(5|'^/ipsmt,ps 



Pt - "^?,ps + ^"^t.PsPt - ^SX^l^Psmt,Ps 



(39) 



23 



to the NLO amplitude and a similar one for the internal antitop propagator. 

The result for the conversion from pole scheme to PS scheme can be equivalently ob- 
tained considering how the top-quark propagator is renormalized and resummed in the two 
schemes. We introduce the short-cut Dt = — ml and -Dt,ps = Vt ~ ''^tps indicate 
the renormalized inverse propagator in the pole and PS scheme (without self-energy re- 
summation). In the pole scheme, the contribution of the hard part of the one-loop QCD 
self-energy to the internal top-quark line can be written as [11] 



2im?(|/t + mt)_ im _ ij^t + m) 
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while the correction to the top-quark line from the mass-renormalization counterterm 6mt 
reads 



6mt 
mt 



(42) 



In the pole scheme one has Smt/rrit = — 5QCD(^t), which leads to the renormalized one-loop 
propagator 



Including in addition the hard part of the electroweak self energy one obtains 

ii^t + mt) 
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In the PS scheme the self-energy correction to the propagator has the same functional 
form as in the pole scheme, but it depends on the PS mass mt^ps 
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Analogously, the mass renormalization gives the contribution 
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The mass counterterm in the PS scheme can be obtained using the relation between the 
pole and PS mass, Eq. (35), and the results given in Eqs. (3.20) and (3.24) of Ref. [11] 
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Thus the one-loop renormalized propagator in the PS scheme reads 
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Remembering that /ips ~ CKsirit and keeping only the parametrically relevant terms one 
obtains 



i(j/t + mt,ps) 
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(48) 



The terms proportional to 5qcb and ^2 ^ scale as ~ cts ~ V6 compared to the leading 
propagator, and can be included perturbatively in the calculation of the cross section. In 

contrast, the term proportional to is a correction of order one, since a s f^'PSf^t, PS / Dt, ps ~ 
1. It thus has to be resummed into the leading-order propagator. Including the contribution 
of the one-loop electroweak self energy, the resummed renormalized propagator reads 



i{i/t + mt, Ps) 



At,ps - ^5f^fjLpsmt,Ps 
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(49) 



with pg = p1 — mjpQ + irrit^ps^t- This confirms that the conversion from pole to PS 
scheme amounts to the replacement At — > A^^ps — ^5PA*ps?7it,ps in the denominator of 
the resummed top-quark propagator and the addition of the term 

l^ps'mt,ps 



"•^pole-i-PS 



Af,ps - '^Sf^fJ'Psmt,PS 
to the NLO amplitude for each intermediate top or antitop quark line present at LO. 



(50) 



3.6 Validity of EFT results 

We conclude this section with a remark on the validity of the EFT results for tt produc- 
tion. As outlined in Section 2, in our effective-theory framework the unstable particle is 
correctly treated by HQET (plus finite- width effects). Strictly speaking, for production of 
a pair of unstable particles this is true as long as the two velocities v, v are relativistic, 
which is the case if ^ 2mt, namely when the partonic centre-of-mass energy \/§ is 
much larger than the pair-production threshold. When \^ — 2mt ~ Tt the two particles 
become non-relativistic, v v (1)0)) ^-^d the correct effective theory to describe the 
pair-production process is non-relativistic QCD (NRQCD). In the expansion by regions 
this is seen as an additional potential region which develops when the heavy particles 
become non-relativistic, and describes Coulomb interactions of the two slow-moving parti- 
cles. This breakdown of the expansion in S was already noticed in the context of the pole 
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approximation for VT-pair production, with the double-pole approximation quickly loosing 
accuracy when approaches 2Miy [70]. While this is relevant for an e~e'^ collider, where 
the centre-of-mass energy can be tuned to probe the threshold region, at hadron colliders 
the physical hadronic cross section is obtained by convoluting the partonic cross section 
and parton luminosities over the full range a/S = [0, ^/s]. For a multi-TeV hadronic c.o.m. 
energy ^/s, and particle masses in the range of 100-200 GeV, the contribution of the thresh- 
old region is expected to be small, so that the error introduced by using the wrong EFT 
description near threshold is negligible. We have checked this explicitly for the total NLO 
cross section in the stable-top approximation and found that the threshold region, defined 
by f < 0.1-0.2, contributes 1-5% of the total hadronic cross section. The Coulomb terms, 
which are the ones incorrectly reproduced by HQET in the threshold limit, are in fact less 
than 1%. 

Below threshold, doubly-resonant configurations do not dominate the cross section 
because the phase space for production of two resonant tops is strongly suppressed for 
\/§ < 2mt. In this region NRQCD should be replaced by a new effective theory describ- 
ing singly-resonant production of a ti* {t*t) pair, which would require the calculation of 
non-trivial QCD corrections to the processes qq — )■ tW~b{iW^b). For this reason, below 
threshold we switch to a Born-level prediction of the top-pair production cross section, and 
set virtual and real corrections to zero. Once again, one expects the error in introducing 
this approximation to be below the accuracy of a few percent which we pursue here. Ac- 
tually, in the results presented next, we impose a physical cut on the final state so that we 
never run into the complications near threshold. 

4 Results 

In this section we present explicit results for a sample of representative kinematical distri- 
butions relevant for top quark pair production studies. Our main purpose here is to discuss 
the importance of off-shell and factorizable corrections by making a direct comparison to 
the NWA approach. The distributions presented in Sections 4.1 - 4.3 are computed in 
the pole-mass renormalization scheme whilst a comparison to results obtained in the PS 
scheme is given in Section 4.4. 

The input values for the SM parameters are given in Table 1. For the convolution of 
the partonic cross section with the PDFs we use the MSTW2008 NLO set [71] for both LO 
and NLO cross sections, taking the corresponding value for for consistency. In addition, 
for all contributions involving the top decay width we use the NLO width, Ff^*"*. Using 
NLO PDFs and r^^° throughout is a choice that has been made to ease identification of 
off-shell effects, ensuring that corrections arising from changes in PDFs or decay widths 
do not cloud the effects introduced by the off-shellness of the tops. Our default choice of 
factorization and renormalization scales is fip = fiji = rrit. Where present, bands around 
LO and NLO off-shell results are obtained by simultaneously varying the renormalization 
and factorization scales in the interval [mt/2,2mj]. 

Throughout the following discussion we assume that the final state contains a b and a 
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Collider 


: Tevatron, ^/s = 1.96 TeV 


PriJ-b, 


) > 15 GeV 
) > 15 GeV 


Pt{1-^) > 15 GeV $T > 20 GeV 
prin > 15 GeV Djct = 0.7 


Mt = 
Mw = 


172.9 GeV 
= 80.4 GeV 


pNLO ^ 1 3gg2 GeV Mz = 91.2 GeV 
Tw = 2.14 GeV ae„ = 0.03394 



Table 1: Example process definition at the Tevatron collider and parameter setup. 

b jet and furthermore that the VT-bosons are perfectly reconstructed. Despite the latter 
not being feasible in the di-lepton channel it allows one to identify features arising from 
the off-shellness of the top quarks that will also be relevant in an experimentally rigorous 
analysis. It is emphasised however, that wherever we discuss an observable related to the 
(anti)top, it is understood this is the reconstructed (anti)top quark, defined in terms of the 
iy+*^~^-boson and (6)6-jet, J^iy 

In order that the EFT approach and counting remains valid the following conditions on 
the reconstructed top-quark momenta, p(t) = p{Jb) + p{W~^) and p(t) = p{Ji) +p{W~), 
are also imposed: 

140 GeV < Minv(t) = \/{p{t)f < 200 GeV, 
140 GeV < Minv(t ) = \l{p{t)f < 200 GeV, 

MiUtt) = \l {p{t) + p{t)f > 350 GeV. (51) 

The experimental setup we consider, typical for a study of top-pair production, is detailed 
in Table 1 and defines the process 

pp JbJi$Tl^ I' + X. (52) 

We note that the jets have been clustered using the standard /cj-algorithm with jet resolu- 
tion parameter Djet = 0.7. 

The implementation of the factorizable corrections to production and decay subpro- 
cesses has been checked in the limit of on-shell top quarks against the di-lepton channel 
implementation in MCFM [63]. Very good agreement was found both for inclusive cross 
sections and for all distributions we have checked. The analytical expressions for the 
virtual non-factorizable corrections have been put through extensive numerical checks and 
corresponding real corrections display the correct independence on the soft FKS parameter 

^cut- 



27 



4.1 Invariant and Transverse Masses 



We start by presenting results for the invariant and transverse masses of the reconstructed 
tops, displayed in Figures 12 and 13 respectively. These observables are expected to be 
sensitive to the off-shellness of the top quarks. The off-shell LO (green band) and NLO 
(red band) results are displayed in the upper plots, as is the NLO prediction in the NWA 
(blue, solid). In order to straightforwardly assess the importance of off-shell effects, non- 
factorizable corrections and sub-leading contributions, we study the quantities 



^NLO, off-s. _ ^NLO, on-s. 
ijNLO, on-s. 



jjNLO, on-s. 



^sub-lead. 
jjNLO,on-s. 



(red, solid), (53) 
(blue, solid), (54) 
(green, dashed), (55) 



in the lower panel of each figure. The non-factorizable corrections have been evaluated at 
the scale yUgoft — frit S, as motivated at the end of Section 3.3. For the case at hand this 
amounts to multiplying o"^^(^ = nit) by the factor Q;s(/isoft)/tts(^t) — 2.24, when 6 = 0.02. 
It would also be desirable to include a resummation of logarithms of /igoft/^t by running 
the hard matching coefficients down to the soft scale using RGEs, however this will not be 
discussed here. 

Firstly we focus on the invariant mass distribution of the reconstructed antitop quark, 
shown in Figure 12.^ An important feature present in the upper plot is that the shapes 
of the NLO off-shell and on-shell distributions are significantly different, particularly near 
and beyond the peak position. This is made clear in the lower inset where one observes 
that off-shell effects are indeed sizeable, exceeding 60% around the peak. However, one 
also notices that they change sign in this region, resulting in the net off-shell effects being 
small due to an averaging out that occurs for inclusive enough observables like the total 
cross section. The lower inset in Figure 12 also shows the effect of strict non-factorizable 
effects and sub-leading contributions, both normalized to the NLO on-shell results. Non- 
factorizable corrections, given by the sum of virtual and real soft contributions, are very 
small, ~ 1% over the whole invariant-mass range considered, except around the pole, where 
they grow to roughly 6-7%. The contribution of sub-leading diagrams is also very small 
close to the peak, where 6^1 and the effective-theory counting is satisfied, but becomes 
sizeable in the tail region for Minv(t) > m^, as one would expect. 

We now consider Figure 13 which displays the transverse mass distribution of the 
reconstructed top quark, defined via 

Mrit) = {{\MJb)\ + \Mn\ + \ET{^i)\f - iMJb)+Mn +M^i)?Y^' (56) 
where Pt(^) is the transverse momentum of the final state k. 



^The LO on-shell result is a delta function centred at rrit and it is omitted in the plot. 
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Figure 12: Invariant mass distribution for the reconstructed antitop quark. Upper plot: 
results for the off-shell LO result (green), off-shell NLO result (red) and on-shell NLO pre- 
diction (blue). Lower inset: relative contributions of total off-shell effects, non-factorizable 
corrections and singly-resonant diagrams. See the text for explanation. 

Once more, it is clear that there are substantial differences between the shape of the 
distribution in the on-shell and off-shell approaches. In particular, the relatively sharp 
edge in the on-shell calculation becomes much less steep once the on-shell assumption is 
relaxed. The size of off-shell effects are quantified in the lower inset where one observes that 
these constitute a negative 2-3% correction to the NWA result, except near and beyond 
the edge (at Mt ~ mt) where much larger effects are present. In this region, the averaging 
effect mentioned earlier is less effective leading to enhanced corrections. The pure non- 
factorizable corrections are negligible for almost the entire range, but do grow to 5-6% 
near the edge. Finally it is clear from the figure that the sub-leading terms are also small 
(~ 2%) for values of Mt below uit and become more important for values beyond the edge. 

The common characteristics shared by the contributions (53) - (55) for this observable 
can be understood by the fact that the distribution in the on-shell approach has a fixed 
edge at mt at LO, whilst at NLO it can only receive contributions beyond nit from events 
involving additional radiation that is clustered into Jb- In the off-shell computation (at both 
LO and NLO) this region receives additional contributions from events where Minv(^) > frtt 
and thus differs from the on-shell result distinctly. It should be noted that in this region 
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Figure 13: Transverse mass distribution for the reconstructed top quark. Upper plot: 
results for the off-shell LO result (green), off-shell NLO result (red) and on-shell NLO pre- 
diction (blue). Lower inset: relative contributions of total off-shell effects, non-factorizable 
corrections and singly-resonant diagrams. See the text for explanation. 

soft-gluon resummation or parton shower effects are also expected to play an important 
role. 

4.2 Individual NLO contributions 

We now examine the relative importance of the different contributions to the cross section 
for the distribution Hx^Jb, Ji) = VT^Jb) + PriJi)- The findings, illustrated in Figure 14, 
are typical of observables that are inclusive in the invariant masses of the top quarks. 

It is clear from the figure that the corrections to the top and antitop decay subprocesses 
are important for the correct normalization of the distribution, each separately correcting 
the LO results by about —10%. The negative sign of the latter means that there is a 
partial cancellation of NLO corrections between the production and decay subprocesses. 
It is evident from the almost flat shape of the decay corrections in the lower inset of 
Figure 14 that these do not noticeably alter the shape of the distribution. Rather, any 
significant correction to the LO shape comes about through the factorizable corrections to 
the production subprocess. The non-factorizable corrections have not been included in the 
plot as they are tiny, as expected by the almost complete cancellation between real and 
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Figure 14: HT{Jb, Ji) distribution. Upper plot: off-shell LO (green, solid) and NLO (red, 
solid) results. The factorizable corrections to the production (blue, solid), top decay (pink, 
solid) and antitop decay (dark green, dashed) subprocesses are also displayed along with 
the tree- level sub-leading contributions (purple, dashed). The ratios of the latter four 
contributions with respect to the LO results are shown in the lower inset (same colours 
and line-styles). See the text for explanation. 



virtual corrections for inclusive observables. In fact, the off-shell effects (also not displayed) 
are , as a whole, also small over the full range of Ht considered, in agreement with the a 
priori expectation that these effects are of order ~ Tt/rrit. The sub-leading contributions 
constitute a 1-2% correction to the LO result (except at very low Ht), indicating that the 
EFT power-counting is working well. 



4.3 Forward-Backward Asymmetry 

Figure 15 displays the differential distribution relevant for the study of the forward- 
backward asymmetry of the top quark, 

da da 

top asymmetry: ^ - (57) 

The important feature to pick out is that off-shell effects amount to a small negative 
correction to the on-shell asymmetry of less than 4% in magnitude. This is fully expected 



31 




-0.001 

-0.002 

-0.003 
0.1 
0.05 




^ 

off-shell correction / on-shell asymm 
sub-leading asymm / on-shell asymm 



0.1 

0.05 



-0.05 



-2 



-1.5 



-0.5 



0.5 



1.5 



Figure 15: Differential forward-backward asymmetry for the reconstructed top quarks, see 
Eq. (57). Upper plot: off-shell (red, solid), on-shell (blue, solid) and sub-leading (purple, 
dashed) asymmetries. Lower inset: ratio of the off-shell correction (red, solid) and sub- 
leading asymmetries (blue, dashed) with respect to the on-shell asymmetry. See text for 
further explanation. 

for an observable that is inclusive over the invariant mass of the top. The sub-leading 
terms display a small asymmetry, giving a ~ 7% correction to the on-shell asymmetry. 

4.4 Pole mass versus PS mass 

In this subsection we examine the effects of using the PS-mass scheme as an alternative 
to the pole scheme. As discussed previously, the PS mass has the advantage of being free 
of non-perturbative ambiguities whilst still being a suitable mass for the effective-theory 
setup as long as the choice /xps ~ cts^t is made. 

We illustrate the above statements in Figures 16 and 17 where the reconstructed top- 
quark invariant mass and pseudo-rapidity are shown respectively. Distributions in the pole 
scheme and PS-scheme with ^pg e {10, 20, 30, 50} GeV at LO (LHS panels) and NLO (RHS 
panels) are plotted. In detail, we fix the numerical value of mt,ps(/ips) to the pole mass 
rrit = 172.9 GeV at next-to-leading order, i.e. including the term ~ in Eq. (35) and 
obtain mt,ps(10) = 172.44 GeV, mt,ps(20) = 171.97 GeV, mt,ps(30) = 171.53 GeV and 

pg(50) = 170.6 GeV respectively. The width of the top is also adapted in accordance 
with rrit at next-to- leading order in as- 

According to our discussion in Sections 2.5 and 3.5, we would expect the suitability 
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Figure 16: Invariant-mass distribution for reconstructed top quark at LO (left panel) and 
NLO (right panel) for pole and various PS-schemes. See the text for explanation. 

of /xps = 50 GeV to be highly questionable, whereas /ips = 10 GeV should be perfectly 
acceptable. To confirm this, firstly, we point out that at LO the distributions in all schemes 
agree very well with each other (within ~ 5% for all bins), even for 'bad' scheme choices. 
This is clear since terms affecting the perturbative stability of a particular mass scheme 
only enter at orders beyond LO. In the RH panel of Figure 16, where the NLO invariant 
mass distributions are plotted, it is evident that the choice /zps ^ 30 GeV leads to serious 
deviations in shape from the pole-scheme curve and to large NLO corrections. This shift is 
due to the term given in Eq. (39) and, to a lesser extent, due to non-factorizable corrections. 
The size of the NLO corrections implies that a NLO description of the lineshape for /xpg > 
30 GeV is not trustworthy. In the same plot however, we observe that sensible scheme 
choices (/ips = 10, 20 GeV) give lineshapes that have a stable perturbative expansion. 
Thus it is perfectly legitimate to use such schemes rather than the pole-mass scheme for 
the extraction of the top mass from a measurement of the invariant mass of its decay 
products. Moreover, at NLO it is only observables that are more exclusive in the invariant 
mass that will display significant differences in shape for different schemes. The latter 
feature is illustrated well for the top pseudo-rapidity in the RH panel of Figure 17, where 
all schemes display the same shape and the NLO corrections are under control. 

In order to investigate the impact of mass-scheme ambiguities on the determination 
of the top mass, we now consider in a toy analysis the invariant mass distribution for 
f^ps ^ {0, 10, 20} GeV, where /ips = of course corresponds to the pole scheme. We want 
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Figure 17: Pseudo-rapidity distribution for reconstructed top quark at LO (left panel) and 
NLO (right panel) for pole and various PS-schemes. See the text for explanation. 



to compare the extraction of the top mass in these schemes at LO and NLO. We stress 
that we only consider effects from scheme changes and ignore all other effects, such as 
colour-reconnection effects and many more. 

Starting at LO, we extract the mass in all schemes by adjusting its value to obtain 
optimal agreement with the measured distribution. Let us assume in the pole scheme the 
measured distribution is matched perfectly for rrit = 172.9 GeV. After the extraction of this 
value for the pole mass we can now convert the latter to the MS scheme to obtain fn = 162.2, 
where the conversion is done at three loop [72] with a crude Pade estimate^ for higher-order 
corrections. However, as we have argued above, we are also entitled to use the PS scheme 
with /ips = 10 GeV or /ips = 20 GeV for such an analysis. Given the results shown in the 
left panel of Figure 16, the extracted values for the masses in these schemes are extremely 
close to the NLO converted masses given previously in this subsection and listed in the 
first column of Table 2. If we now convert these values to the MS scheme with precisely the 
same procedure used for the pole scheme, we obtain the values for fn listed in the second 
column of Table 2. These values differ by 600 — 800 MeV from m = 162.2. Alternatively, 
if we convert the extracted values of mt_ps(10) = 172.44 GeV and mt^ps(20) = 171.97 GeV 
back to the pole scheme, we obtain the values rrit = 173.5 GeV and rrit = 173.8 GeV 

^We emphasise that the precise method employed to estimate effects of higher order corrections does 
not play a major role here. For all numbers quoted in this subsection, including Table 2, the estimate of 
the error in the conversion is < 100 MeV. 
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Figure 18: Invariant-mass distribution for reconstructed top quark at NLO for pole and 
various PS-schemes. See the text for explanation. 

respectively, as listed in the third column of Table 2. This conversion is done at 0{aj), 
again supplemented by a crude Pade approximation to estimate higher-order effects. We 
stress that all values given in the left panel of Table 2 are of course consistent at 0{as) 
and differ only by terms that are formally of higher order. The crucial point is that the 
numerical effect of these higher-order terms lead to an ambiguity of 500 — 900 MeV in the 
LO extraction of the top mass from the invariant mass distribution. 

We can now repeat this exercise at NLO where we begin by stressing that NLO in this 
context implies the inclusion of NLO corrections to the propagation of the top quark, in 
addition to the factorizable and non-factorizable corrections. Simply taking into account 
NLO corrections to the production subprocess alone does not improve the LO toy analysis 
presented above. We again assume that the experimental measurement of the invariant 
mass is perfectly matched in the pole scheme by setting rrit = 172.9 GeV. To extract the 
PS-mass from this measurement we would have to perform a best-fit analysis in the PS- 
scheme and extract the best value for the top mass. We have found that this best value 
is very close (but not exactly identical) to the value obtained from two-loop conversion of 
the pole mass to the PS mass, listed in the first column in the right panel of Table 2. In 
Figure 18 we show the corresponding distributions to confirm that they are indeed very 
close. 

As in the LO case we then convert the extracted values for the mass to the MS scheme 
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LO 






NLO 




A*ps 




m 


mt 


m^xp 


m 


mt 





172.9 


162.2 


172.9 


172.9 


162.2 


172.9 


10 


172.4 


162.7 


173.5 


172.2 


162.4 


173.3 


20 


172.0 


163.0 


173.8 


171.5 


162.5 


173.4 



Table 2: Extraction of the top mass in various schemes at LO (left panel) and NLO (right 
panel). All numbers are in units of GeV. 

(second column) and back to the pole scheme (third column). Not surprisingly, the ambi- 
guities in the extraction of the mass has decreased compared to the LO case and is now 
about 200 — 300 MeV for the determination of m and 400 — 500 MeV for the determination 
of the pole mass. 

Let us conclude this subsection with a few comments regarding this toy analysis. First 
and foremost, this is of course only a basic investigation. Apart from taking into account all 
partonic channels, a fully rigorous analysis would have to include numerous effects beyond 
fixed order in perturbation theory. However, the main result that there is an existing 
additional ambiguity due to the scheme dependence, is unlikely to change drastically. In 
particular, it is not justified to blindly identify the extracted mass of the top quark as 
the pole mass. The size of this additional error will depend most of all on whether the 
analysis is done at LO or NLO. It also depends on what values of /xps are considered 
acceptable - increasing the range in /xps increases the spread in the extracted mass. We 
find that the ambiguity in the (indirect) determination of fn is smaller than that present 
in the determination of the pole mass. This could be due to the better convergence of the 
associated perturbative series relating the different masses, which in turn is related to the 
renormalon ambiguities inherent to the pole mass. 

5 Conclusions 

We have presented an approach that allows the computation of cross sections involving 
heavy unstable particles, formalising the physical picture of production and subsequent 
decay. Our approach makes use of techniques common to effective-theory calculations, 
organizing the amplitude into contributions from matching coefficients and from dynamic 
soft degrees of freedom in the effective theory. This spht allows one to disentangle and 
separately study effects from widely different scales. Compared to computations using 
the complex-mass scheme, our approach leads to calculations that are considerably sim- 
pler technically. However, our results do have the disadvantage of being valid only in the 
resonant region. To relax this constraint within the effective theory approach, we would 
have to match our results to a strict fixed-order calculation (without any self-energy re- 
summation) outside the resonance region. Similarly, if we wanted to extend the results to 
the threshold region, we would have to match our effective theory to another NRQCD-like 
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effective theory that describes a top pair near resonance. 

The approach apphed in this paper to top pair production from initial state quarks is 
a generahzation and development of an approach applied earlier to single-top production. 
In particular, the real corrections are expanded at the integrand level so that no further 
expansions after phase space integration are required. 

Generically, off-shell effects are found to be small and only have an impact near kine- 
matic edges where distributions are not inclusive in the invariant mass and, therefore, the 
corresponding cancellations are not complete. As a toy application we have investigated 
the impact of mass-scheme ambiguities for the extraction of the top mass from the invari- 
ant mass distribution. In a simple analysis we found scheme ambiguities of 500 — 900 MeV 
and 400 — 500 MeV for the determination of the pole mass using LO and NLO calculations 
respectively. Future apphcations of the EFT method include the possibility of studying 
resummation of logarithms of ratios of widely separated scales for fully differential cross sec- 
tions. The fact that effective theories disentangle effects from different scales is a promising 
starting point for such an investigation. 
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